Timescale for equilibration of N/Z gradients in dinuclear systems 
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Equilibration of N/Z in binary breakup of an excited and transiently deformed projectile- like 
fragment (PLF*), produced in peripheral collisions of 64 Zn + 27 Al, 64 Zn, 209 Bi at E/A = 45 MeV, 
is examined. The composition of emitted light fragments (3<Z<6) changes with the decay angle of 
the PLF* . The most neutron-rich fragments observed are associated with a small rotation angle. A 
clear target dependence is observed with the largest initial N/Z correlated with the heavy, neutron- 
rich target. Using the rotation angle as a clock, we deduce that N/Z equilibration persists for times 
as long as 3-4 zs (lzs = 1 x 10 _21 s = 300 fm/c). The rate of N/Z equilibration is found to depend 
on the initial neutron gradient within the PLF*. 
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The broad impact of the density dependence of the nu- 
clear symmetry energy makes it a topic of considerable 
interest. Whether the asymmetry term in the nuclear 
equation of state follows a stiff or a soft dependence on 
density determines the composition of a neutron star's 
crust, and the conditions under which a supernova explo- 
sion occurs [IHH . In the case of nuclei, the stability of the 
heaviest elements 0], and the existence of neutron skins 
H, also depend on this quantity. One means of investi- 
gating the density dependence of the symmetry energy is 
by measuring the N/Z equilibration in a dinuclear system 
[Q-Q. Although past studies have principally focused on 
the N/Z equilibration that occurs between the projec- 
tile and target nuclei in a collision, such an approach is 
fundamentally limited by the short contact time between 
the two collision partners. This contact time, which is 
approximately 100 fm/c at intermediate energies, inher- 
ently limits the degree of equilibration that can be at- 
tained. 

Another opportunity to investigate N/Z equilibra- 
tion is the dynamical binary breakup of an excited 
and transiently deformed nucleus produced in the semi- 
peripheral collision of two heavy-ions at intermediate en- 
ergies (E/A=20-100 MeV) %. At these energies, the 
collision of two heavy-ions at peripheral and mid-central 
impact parameters leads to the exchange of charge, mass, 
and energy between the projectile and target nuclei. 
Prior experimental and theoretical work demonstrates 
that several mechanisms can contribute to the ternary 
breakup of the transient system into two heavy remnants 
(projectile and target) together with one intermediate 
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mass fragment For many such breakups, the ex- 

istence of a short-lived neck joining the projectile and 
target plays a crucial role. As the two heavy remnants 
separate the neck ruptures leaving it still attached to one 
of the two remnants. The binary system consisting of the 
neck and associated remnant continues to undergo neu- 
tron and proton exchange as it rotates. If the neck is 
associated with the projectile remnant, we designate the 
joint neck-heavy remnant system as the projectile-like 
fragment (PLF*). In the case of ternary decays a second 
rupture of the neck occurs producing the intermediate 
mass fragment. Statistical decay from the equilibrated 
projectile and target remnants also occurs. 

Equilibration of N/Z within the PLF* has been ob- 
served to persist for timescales as long as t w 3zs (lzs = 
1 x 10 _21 s) [l3|. In the present work we further charac- 
terize this approach, extending it to a significantly lighter 
projectile. We explore the impact of the N/Z of the tar- 
get on the composition of the fragments produced and 
their subsequent evolution. We also demonstrate that 
this equilibration process exists when the projectile and 
target nuclei are the same nucleus. 

The experiment was conducted at the Cyclotron In- 
stitute at Texas A&M University, where a beam of 64 Zn 
ions was accelerated to E/A = 45 MeV with an aver- 
age beam intensity of «2xl0 8 p/s. The beam impinged 
on 27 Al, 64 Zn and 209 Bi targets with thicknesses of 13.4, 
5, and 1 mg/cm 2 respectively. Although the experimen- 
tal details have been previously published [l4[, they are 
summarized below for completeness. The array FIRST 
[la] .which subtended the angular range 4.5° < < 
27° was used to identify charged products produced in 
the reaction. In the angular range 4.5° < 8i a b < 7°, the 
forward telescope in FIRST provided identification by 
atomic number of all products up to Z=30 and isotopic 
information for Z < 12. The second telescope in FIRST 
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FIG. 1: (Color online) The angular distribution of binary 
splits (Zl-Zh) for Zl=4, representative of other fragments, 
is shown. Data for the 64 Zn, 27 A1, and 209 Bi targets are rep- 
resented by the black line, blue open symbol and red closed 
symbol histograms respectively. 



, which subtended the angular range 7° < 9i ab < 14°, 
provided Z identification for Z < 22 and A identification 
for Z < 8. The third telescope in FIRST (14° < 9 tab < 
27°) provided Z identification for Z < 12 and A iden- 
tification for Z < 7. The high segmentation of FIRST 
provided an angular resolution of ±0.05°(4.5° < 9i ab < 
7°), ±0.44°(7° < 9 lab < 14°) and ±0.81°(14° < 9 lab < 
27°) in polar angle and ±11.25° in azimuthal angle. The 
energy resolution obtained was approximately 1%. 

In order to focus on binary decays, events were selected 
in which two fragments (Z > 3) were detected within the 
laboratory angular range 4.5° < 9i ab < 27°. These two 
fragments were distinguished from each other by their 
atomic number, with the larger (smaller) atomic frag- 
ment designated as Zh (Zl). We ensured that the PLF* 
under investigation comprised a large fraction of the ini- 
tial projectile by requiring that the events selected had 
Zh > 11. Events selected in this manner corresponded 
to approximately 14% of the measured yield in which one 
fragment with Z>11 was detected. 

It has previously been established that an instructive 
quantity for the binary decay of the PLF* is the angle 
between the direction of the two fragments center-of-mass 
velocity, v c ,m , and their relative velocity, vrel, denned 
as vh-vl [UllH]- We construct the angle a, with 

COs(a) = V c .ra. ■ V RE L / (\\v c .m. \\ \\VREL ||) 

as indicated within the inset of Fig. Q] Consequently, 
aligned decays with Z^ emitted backward (forward) of 
Zh corresponds to cos(a) — 1 (-1). Momentum correla- 
tions observed between Zh and Z l reveal that these two 
fragments originate from a common parent as evident in 
Fig. 2 of [17| . This parent nuclear system comprised of 



Zh and Zl is designated as the PLF*. 

The angular distribution for Z^=4 fragments is pre- 
sented in Fig.[T]for the 64 Zn, 27 A1, and 209 Bi targets used. 
Consistent with previous work, all the angular distribu- 
tions manifest a peak at backward angles cos(a)>0.5. 
This preferential backward decay of the PLF* is well es- 
tablished, and has been interpreted as the aligned dy- 
namical decay of the excited and transiently deformed 
PLF* [H [H, [l]|. The backward peaking of the angular 
distribution can be understood as the dynamical binary 
splitting of the PLF* on a timescale that is short rela- 
tive to the rotational period of the PLF*. To compare 
the shape for the different targets, all three distributions 
have been normalized in the interval -l<cos(a)<0. We 
chose this region for normalization since it corresponds 
to forward statistical emission from the PLF*. This for- 
ward statistical decay is long-lived relative to backward 
emission and hence is less coupled to any dynamics re- 
sponsible for the formation of the PLF* . For each target, 
the angular distributions exhibit the same shape, mani- 
festing a distinct preference for aligned decay of the PLF* 
with the Zl fragment oriented towards the target. The 
shape of the distribution for cos(a)<0 provides an indi- 
cation that the angular momentum of the decaying PLF* 
is relatively small. In contrast to previously studied sys- 
tems |l3[, the yield does not increase near cos(a)=-l. 

Given the normalization at forward angles, the similar- 
ity of the distributions for the three targets at backward 
angles, cos(ev)>0, is striking. This similarity suggests 
that while the probability of forming the elongated and 
excited PLF*, as well as its composition, may depend 
on the target, the relative probability of its subsequent 
decay is essentially independent of the target. 

We next examine whether the composition of the Zl 
fragment changes as a function of rotation angle. In 
our initial work, which analyzed the reaction 124 Xe + 
112 < 124 Sn, we observed that the (N)/Z of the Z L frag- 
ment decreased as the Z l-^h system rotated [l3| . As the 
Zn-like PLF* in the present work is considerably smaller 
than the Xe-like PLF* , it was unclear whether the behav- 
ior previously observed would also exist for the smaller 
PLF*. Depicted in Fig. His the (N)/Z of the Z L fragment 
for Be (panel a), B (panel b), and C (panel c) fragments. 
For each Z l shown the impact of the three different tar- 
gets is also presented. The average (N)/Z for a given 
Zl is deduced by averaging the neutron number for the 
different isotopes measured. A common feature of all 
the data is that the largest value of (N)/Z is associated 
with cos(a)=l, namely backward emission. As the Z^- 
Zh system rotates, (N)/Z of the Z^ fragment decreases 
corresponding to a net loss of neutrons by the Z^ frag- 
ment. In the case of the Be fragments, this dependence of 
(N)/Z on cos(a) is clearly apparent even for the lightest 
target, Al. For all three fragments shown the magnitude 
of (N)/Z is largest for the Bi target. We attribute this 
large value of (N)/Z for cos(a)~l in the case of the Bi 
target to the preferential pickup of neutrons by the PLF* 
from the Bi target with its N/Z=1.51. In contrast, the 
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FIG. 2: Average neutron to proton ratio for selected Zi as 
a function of the decay angle. The ratio for the 64 Zn, 209 Bi, 
and 27 Al targets is represented by the closed circle, open circle 
and open triangle respectively. The dashed line for each Zl 
represents the N/Z of the 64 Zn beam. 



64 Zn and 27 A1 targets with N/Z = 1.13 and 1.07 do not 
present a neutron-rich reservoir from which the PLF* can 
pick up neutrons. 

Shown in the right hand scale of Fig. [2] is the (N) of 
the Zl fragment. In the case of the Be fragments, for the 
Bi target (N) decreases from 5.2 to 4.05 a net change of 
over one neutron. For the Zn and Al targets, a somewhat 
smaller net decrease of 0.6 - 0.7 in neutron number is ob- 
served. For the Bi target, the change in (N) for Z^=5 
and Zl=6 is «0.2. The change in (N) for Li fragments 
and the Bi target (not shown) is also «0.2, comparable to 
that of B and C fragments. The larger change observed in 
the case of Be fragments can be qualitatively understood 
as being due to the absence of 8 Be fragments. Since the 
isotopic distribution for all the fragments with Zl=3,5, 
and 6 has a value of (N)/Z>1, it is reasonable to expect 
that this is also the case for Be fragments. The decay 
of 8 Be into two alpha particles removes these fragments 
from the measured isotopic distribution thus artificially 
increases the value of (N) observed for Be at backward 
angles. In effect, the absence of 8 Be acts as an amplifier 



for the change in (N) by emphasizing the importance of 
the extremes of the isotopic distribution. This conclu- 
sion is supported by our re-analysis of carbon isotopes in 
which we eliminate 12 C from the isotopic distribution [20|. 
For this reason, we have elected to present the Be data 
without correcting for the absence of 8 Be. 

The physical picture that emerges is one in which the 
N/Z of the dinuclear PLF* is established through its in- 
teraction with the target. Preferential transfer of neu- 
trons from a neutron-rich target such as Bi results in a 
neutron- rich PLF*. As the nascent Z^ fragment is ori- 
ented towards the target-like fragment, it is the primary 
beneficiary of the transferred neutrons. In addition, even 
for a symmetric projectile-target collision, the density de- 
pendence of the symmetry energy leads to neutron en- 
richment of the low-density neck [14]. The result is an 
initial N/Z gradient within the PLF*. As time passes, 
these additional neutrons in the Z^ fragment are dissi- 
pated. Whether this preferential neutron transport out 
of the Z l fragment occurs into the Z# fragment or into a 
low-density neck region connecting the Z^ and Z# frag- 
ments is presently unclear [jjj]. It should be clear that 
transfer of both neutrons and protons occurs between the 
Zjj and Zi fragments. Our selection of a particular Z^ 
fragment in this analysis precludes us from examining the 
net proton exchange. 

If the decrease in (N) /Z with cos(a) can be understood 
as the preferential transport of neutrons out of the Z^ 
fragment, one might expect that the shorter the contact 
time between the Z^ and Z# fragments the less likely it 
is that the initial (N)/Z is decreased. Dynamical split- 
ting of the dinuclear Z^-Z^ system can be viewed as a 
dynamical fission process in which the reaction dynam- 
ics provides collective motion along the separation axis 
of the Zl-Zh system [HEEE3- 

Within such a picture 
we expect that the shortest times (dynamical ruptures) 
are associated with the largest relative velocities and the 
longest times are associated with smallest (Coulomb bar- 
rier) relative velocities. 

In order to explore the dependence of (N)/Z on both 
cos (a) and vrel, we present the dependence of (N) of 
Be fragments in velocity space in Fig. [3] In this figure, 
the dependence of the average neutron number, (N), of 
Be fragments on the transverse (vrel.±) and parallel 
(vrel,\\) components of vrel is depicted. The parallel 
and transverse components of vrel are calculated with 
respect to the center-of-mass velocity of the Z^-Z# sys- 
tem. For reference, relative velocities between 1.5 and 
5.5 cm/ns are indicated as dotted circles while the an- 
gular cuts over which the average neutron number was 
calculated are represented by dashed lines. 

For the 209 Bi (panel a) target a systematic behavior of 
(N) of the Be fragment is observed. As one rotates clock- 
wise in the two dimensional velocity space, i.e. increas- 
ing rotation angle a, the value of (N) decreases. For the 
largest vrel, the (N) of the Z^ fragment decreases from 
5.5 to 4.46, a change of » 1 neutron as the Zh-Zl system 
rotates by a quarter turn. For the two most backward an- 
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FIG. 3: Dependence of the average neutron number on the 
parallel and transverse components of vrel for Z_l = 4 for 
the reactions 64 Zn + 209 Bi, 64 Zn.(See text for details) 



gle bins one observes that (N) decreases with decreasing 
y rel,\\ (from 5.5 to 4.6), while for larger rotation angles, 
(N) increases with decreasing vrel,\\ (from 4.46 to 4.58). 
From this trend, one would predict that (N) and hence 
(N)/Z would decrease with increasing vrel for forward 
emission. This expectation is confirmed for this system in 
agreement with previous observation (l3| . These trends 
are also observed for the 64 Zn (panel b) target although 
the magnitude of the change is slightly smaller than in 
the 209 Bi case. 

To extract the time dependence of the (N) of the Zj, 
fragment, we utilize the rotation angle of the Z^-Zr din- 
uclear system as a clock such that the rotation time, t 
is given by: t — ajuj where u> is the angular frequency. 
Hence, the quantity to be determined is the angular fre- 
quency which is given by: uj — (Jft)/I e // where J is 
the angular momentum and I e / / is the moment of iner- 
tia for the dinuclear system. The angular momentum of 
the dinuclear complex is determined by utilizing a sim- 
ple model that describes the statistical decay of a ro- 
tating system, appropriate for forward emission. This 
model was chosen because it provides a simple way to 



include the effects of thermal excitation and collective 
rotation but may have some limitations. The magnitude 
and direction of the velocity of the PLF* are determined 
by sampling the experimental data for forward emission. 
The magnitude of the relative velocity vector between 
Zr and Z^ is taken from the Viola systematics [2l| with 
a width provided from the experimental data. The in- 
plane and out-of-plane components of vrel are calcu- 
lated relative to the plane defined by the PLF* and beam 
direction. The distribution of the out-of-plane emission 
of Z^ is taken as: P(smcj))=Aexp(-k 2 sin 2 (j>) where <j> is 
the out-of-plane angle, k represents the width of the 
distribution, and A is a normalization constant. The 
model predictions have been filtered by the detector ac- 
ceptance and compared to the experimentally measured 
angular distributions presented in Fig. Q] Comparison 
of the measured and predicted distributions for differ- 
ent values of k indicates that the magnitude of k is w 
0.5. Within the framework of a fissioning nucleus, the 
parameter k can be related to angular momentum [22j : 
J 2 = (2k 2 I e f fT) /ft 2 where T is the temperature and I e f f 
is calculated as: I e ff = ^MR 2 Fj. The mass, M, is ap- 
proximated as: M = moc 2 ApLF* where moc 2 is the rest 
mass of the nucleon and A PLF - = (4) ■ , , Z PLF ,. 

r ^ r \ £ / projectile 

The effective radius of the dinuclear configuration is given 
by R 2 F/ with R=r A 1 / 3 and the deviation from a sphere 
accounted for by Fj [23j . The value of the radius constant 
ro is taken as 1.2 fm. As F/ has not been calculated for a 
system as light as the PLF* under consideration, we use 
the published value for the significantly heavier nucleus, 
149 Tb 23]. Assuming a temperature of T=3-5 MeV for 
the system undergoing binary decay, we calculate an an- 
gular momentum J=6±l ft. 

The timescale deduced in this manner is shown in 
Fig. |31 It should be noted that the timescale deduced 
(t < 3 zs) is consistent with previously published results 
[13, [24|,|25|]. For reference, the angular velocity calculated 
for this light dinuclear complex is 0.4-0.5 x 10 21 rad/sec. 

Having associated the rotation angle with a timescale, 
it is now possible to observe two timescales evident in 
Fig. [3] The first observation is that for this system even 
for times as long as 3 zs, i.e. 900 fm/c, the (N) of the Zl 
fragment is still changing indicating that N/Z equilibra- 
tion is a slow process. In addition, operating on a faster 
timescale of w 1 zs, the vrel dependence observed for 
the most backward angles is overpowered by the Coulomb 
effect that characterizes forward emission 26] . The pat- 
tern observed for Z^=4 in Fig. [3] indicates that even for 
backward angles a correlation exists between the rota- 
tion angle dependence and the vrel dependence for (N). 
However, disentangling the intrinsic N/Z gradient from 
the Coulomb contribution at backward angles requires 
knowledge of the detailed configuration of the dinuclear 
system. The trajectory of the Zl fragment depends on 
the motion of the Zl fragment relative to both the target- 
like and Zr fragments which perturbs the intrinsic N/Z 
pattern. This disentanglement is beyond the scope of the 
present work. We therefore examine the dependence of 
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FIG. 4: (Color online) Dependence of the average differential 
neutron number on time for Zl = 4, 5 and 6 for the reactions 
64 Zn + 64 Zn, 209 Bi and 27 Al.(See text for details) 



(N) on time, integrated over vrel- 

Shown in Fig. |4] is the time dependence of (N) for 
Z^=4, 5, and 6 fragments. To compare the change in 
(N) for different targets and different Zl fragments, we 
subtracted the extrapolated value of (N) at t=0. This 
extrapolated value was determined by performing a lin- 
ear fit of (N) versus time. It is interesting that the data 
shown in Fig. |4] clearly fall into two groups. While Z^=4 
exhibits a strong dependence of A(N) on time, Zl=5 and 
6 fragments manifest a much weaker dependence. In the 
case of Z^=4 fragments, a strong and clear target depen- 
dence is evident. Close examination of the data for Zl=5 
and 6 fragments reveals that the same target dependence 
exists though it is smaller in magnitude. For Zl = 4, the 
Bi target is associated with the largest change in (N), «1. 
The Zn and Al targets exhibit smaller changes of 0.75 and 
0.6 respectively. The horizontal error bars shown in the 
figure are primarily governed by an estimated uncertainty 
of AJ=l?i. It should be noted that the magnitude of the 
change in (N) is considerably larger than that previously 
reported for the Xe + Sn system [ill ]. 

For Be fragments the large slope of A(N) versus time 
indicates significant equilibration occurs on the timescale 
of 3 zs. One can imagine that the rate of equilibration 
is governed by the initial difference in N/Z between the 
Zl and Zh fragments, i.e. a gradient in N/Z within the 
PLF* . Confirming this perspective is the behavior of Be 



fragments for different targets. The data for the Bi target 
manifests the largest slope and the one for the Al target 
the smallest slope corresponding directly to the neutron- 
enrichment of the Z l fragment as evident in Fig. [5] This 
target dependence suggests a fundamental result. The 
neutrons transferred from the target do not equilibrate 
within the PLF* prior to it attaining the dinuclear con- 
figuration. If they did, all three targets would manifest 
the same equilibration rate. Thus the experimental data 
directly indicate that the equilibration rate of (N) is gov- 
erned by the initial N/Z gradient in the dinuclear system. 

In summary, to investigate the N/Z equilibration 
timescale for a dinuclear complex, we have focused on the 
binary decay of a PLF* . For this system, we have shown 
that the isotopic composition ((N)/Z) of the fragments 
depends on the rotation angle of the dinuclear complex. 
Moreover, the N/Z of the Zl fragment depends on the 
neutron- richness of the target. While we observe these 
effects for Z/,=4, 5, and 6 fragments, they are strongest 
in the case of Be fragments. This large magnitude for Be 
fragments can be qualitatively understood as due to the 
absence of 8 Be in the isotopic distribution making Be a 
sensitive probe of the equilibration process. Within our 
analysis, the N/Z equilibration is related to the change 
in average neutron number, A(N), experienced by the 
Zl fragment. Using the rotation angle of the dinuclear 
complex we deduce the timescale on which the N/Z equi- 
libration occurs. The equilibration persists for as long 
as 3-4 zs. Clear evidence for a target dependence of 
the equilibration rate is observed in the case of Be frag- 
ments. The largest equilibration rate (Bi target) is asso- 
ciated with the largest initial neutron-enrichment of the 
Zl fragment. These differences in equilibration rate can 
be related to different initial N/Z gradients within the 
dinuclear configuration for the different targets. While 
qualitatively the physical picture associated with these 
observations is clear, a more detailed understanding of 
the N/Z equilibration will require comparison with a the- 
oretical model that describes nucleon transport within 
the dinuclear complex. 
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